Ballistic spin field-effect transistors: Multi-channel effects 
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^ i Abstract 

We study a ballistic spin field-effect transistor (SFET) with special attention to the issue of multi- 
channel effects. The conductance modulation of the SFET as a function of the Rashba spin-orbit 
■ coupling strength is numerically examined for the number of channels ranging from a few to close 

<p : 

to 100. Even with the ideal spin injector and collector, the conductance modulation ratio, defined 
as the ratio between the maximum and minimum conductances, decays rapidly and approaches one 
with the increase of the channel number. It turns out that the decay is considerably faster when the 
Rashba spin-orbit coupling is larger. Effects of the electronic coherence are also examined in the 

o 

multi-channel regime and it is found that the coherent Fabry-Perot-like interference in the multi- 
channel regime gives rise to a nested peak structure. For a nonideal spin injector/collector structure, 

. which consists of a conventional metallic ferromagnet-thin insulator-2DEG heterostructure, the 

OO , 

S^Q , Rashba-coupling-induced conductance modulation is strongly affected by large resonance peaks 

o 

that arise from the electron confinement effect of the insulators. Finally scattering effects are 

\o : 

briefly addressed and it is found that in the weakly diffusive regime, the positions of the resonance 



peaks fluctuate, making the conductance modulation signal sample-dependent. 



*Electronic address: jsjeong@postech.ac.kr 
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I. INTRODUCTION 



At an interface of a semiconductor heterostructure, electrons confined to the interface 
form a two-dimensional electron gas (2DEG). The Rashba spin-orbit (RSO) coupling 
arises in the 2DEG when the interface confinement potential breaks the structural inversion 
symmetry. It is demonstrated that the strength of the RSO coupling can be modulated by 
a gate voltage that affects the asymmetry of the confinement potential. The RSO coupling 
can be a useful tool for spintronic applications in semiconductors [3] . One of representative 
examples is the spin field-effect transistor (SFET) proposed by Datta and Das j^, which 
is based on the spin precession by the RSO coupling within the 2DEG and on the spin 
preparation and detection by two spin-selective electrodes (injector and collector). 

Despite intensive experimental efforts, the SFET has not been realized yet and there are 
numerous reports addressing practical problems for the realization of the SFET such as the 
spin injection/detection [5|,|6| and the spin relaxation 

BBS 

. In this paper, we focus on the 
issue of the multi-channel effects. While a narrow 2DEG with only one transport channel is 

n 

an ideal environment for the SFET operation 10], a single-channel system is rather difficult 
to realize in experiments. For a 2DEG with the Fermi wavelength of 100 A, for instance, the 
width of the system needs to be smaller than 100 A for the system to be in the single channel 
regime. Preparing a 2DEG with its width less than 100 A is, though not impossible, quite 
demanding. The experimental difficulty is further enhanced when it is taken into account 
that for the transistor operation, the system needs to be longer than a certain minimum 
length for the spin precession angle to be of order ir at least. For the 2DEG formed at 
the InGaAs/InAlAs interface, for instance, the minimum length was estimated to be 0.67 
/im 4|. Taking into account the experimental difficulty in the preparation of a sufficiently 
long single-channel system, a multi-channel system is thus a more practical test ground for 
the SFET. 

Behaviors of a multi-channel SFET may be crucially affected by the inter-channel cou- 
pling. In Ref. J], it was argued that effects of the inter-channel coupling will be negligible 
and a multi-channel SFET will behave as a collection of uncoupled single-channel SFETs 
when a dimensionless number 2m*aw/h 2 is sufficiently smaller than 2tt, where w is the 
width of the 2DEG, m* is the effective mass of the electron in the 2DEG, and a is the RSO 
coupling parameter. In this paper, we examine systematically behaviors of a multi-channel 
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;o 100 and with a from 10 12 eVm to 



12 



13, 



141 ] . For the examined range, 



SFET with the number of channels from a few to close 
10 -10 eVm, covering most reported values of a 
the value of 2m*aw/h 2 varies from approximately 0.03 to 50, including the range where the 
inter-channel coupling is not negligible. Our result thus provides a systematic investigation 
of the multi-channel effects. We also examine effects of the electronic coherence in the multi- 
channel regime. For a single-channel SFET, it was recently demonstrated 
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that the electronic coherence makes the conductance deviate from the sinusoidal dependence 
on a. We find that the electronic coherence in the multi-channel regime can give rise to 
a nested peak structure. Our analysis indicates that the nested peak structure is due to 
the coherent Fabry- Perot- like interference in the multi-channel regime. We remark that the 
nested peak structure was not found in previous studies of a multi-channel SFET 

since effects of the electronic coherence were not properly taken into account. 

The paper is organized as follows. Section [TT1 reports the numerical conductance calcula- 
tion results of the conductance of multi-channel SFETs as a function of the RSO coupling 
parameter a. SFETs with various channel numbers are examined. When the spin injector 
and collector in SFETs are ideal, it is found that the conductance modulation ratio between 
the maximum and minimum conductances reduces as the number of channel increases. The 
decaying rate of the ratio is significantly larger for SFETs with larger a. In order to gain 
insights into the numerical calculation results, energy dispersion relations and spin configu- 
rations of energy eigenstates are examined. Section IrTTl focuses on the conductance properties 
of SFETs with rather small number of channels that exhibit a large conductance modulation 
ratio. It is found that the electronic coherence gives rise to a nested peak structure. Effects 
of an in-plane magnetic field are also addressed briefly and it is verified that the magnetic- 
field- induced peak splitting, reported previously for a single-channel SFET jlij, persists in 
the multi-channel regime. Section IIVI addresses briefly behaviors of multi-channel SFETs 
equipped with nonideal spin injectors and collectors, which consist of conventional metallic 
ferromagnet-thin insulator-2DEG hybrid structures. Scattering effects are also addressed 
briefly. Section IVl summarizes the paper. 
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II. CONDUCTANCE OF MULTI-CHANNEL SFET 

Figure ^ shows a schematic drawing of a SFET with finite width w. Within the 2DEG 
in the xy-plane, the effective Hamiltonian of electrons reads 



where p = {p x ,P y ) is the momentum operator within the 2DEG, V c (y) is the transverse 
confinement potential with V c (y) = for < y < w and V c (y) = oo otherwise. Here 
a = (a x , a y , a z ) is the Pauli spin operator, and z a unit vector perpendicular to the 2DEG. 

In order to focus on the issue of the multi-channel transport, we assume situations where 
effects of other practical problems are minimized; scattering by impurities and phonons is 
ignored in if 2 D and ideal spin injection and detection are assumed. To be specific, an electron 
in the injector and collector is assumed to be subject to the following Hamiltonian, 



where B ex is the effective exchange energy in the injector/collector, g is the Lande's g- 
factor, and fiB is the Bohr magneton. For sufficiently large gfisB ex (>0), electrons in 
the injector and collector are 100% spin-polarized along the (+x)-direction. Note that 
with this choice of H- m y co [, the Fermi wavelength of (+x)-polarized electron within the 
injector /collector matches perfectly with that within the 2DEG when a = and thus the 
so-called conductance mismatch problem for the spin injection/detection does not appear 
for the 2DEG-injector(collector) interface. For nonzero a, there is a weak mismatch in the 
Fermi wavelength but for the range of a values examined below (a < 10~ 10 eVm), the 
conductance mismatch problem turns out to be a weak effect (see Appendix El for details). 



FIG. 1: Schematic diagram of a ballistic spin field effect transistor consisting of two-dimensional 
electron gas (2DEG) of finite width w and two spin-selective electrodes (injector and collector), 
whose spins are assumed to be fully polarized along the (+a;)-direction. The strength of the Rashba 
coupling can be tuned by external gates (not shown). 
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A. Numerical conductance calculation 



For the numerical conductance calculation of the mult i- channel SFETs, we use the tight- 
binding (TB) Hamiltonian 0,Ei| H TB = #™ + #™ + #™+#™ pling , where Hj B amounts 
to the TB approximation of the Hamiltonian H 2 t> [Eq. (JU] for the 2DEG, 
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where Ci j Sz is the annihilation operator of an electron at = a(ix + jy) with spin s z 
along the z-axis, a is the lattice spacing used for the TB approximation, and N x and N y are 
related to the length L and width w of the 2DEG via a(N x + 1) = L and a(N y + 1) = w, 
respectively. With the choice t = h 2 /2m*a 2 and A = a/2a, Eq. Q reduces to Eq. ((TJ) in the 
limit a — > 0. The TB Hamiltonians for the injector and collector are similarly given by 

— OO Ny 
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Y Y ( C l+l, j>Sx =lCi,3,s m =l + H.c.) + Y Y ( C h+l,s x =lCi,j,s m =l + H.c.) 



i=N x +l j=l 



i=N x +l j=l 



where c,- 



Hj, Sx =±i = (q JjSz= i ± c i; j >Sz= -i)/\r2. Note that the electron operators for s x = —1 do 
not appear in -ff™ col , since the injector /collector is assumed to be 100% spin polarized along 
(+x)-direction. The coupling Hamiltonian between the 2DEG and the injector/collector is 
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given by 

Ny 

-^coupling = ~t ( C l,i, Sa; =l c 0,j,^=l + C ] Nx+l j Sx=l C Nx j, Sx= i + H.C.J . (6) 
3=1 

Note that electrons whose spin is pointing along the (— x) -direction are not allowed to hop 
between the 2DEG and the injector (collector). Note also that the hopping parameter t in 
iJ™ col is the same as that in if™. Thus when effects of a (or A) on the Fermi wavelength 
is not significant, the conductance mismatch problem should be minimal (see Appendix 1X1 
for detail), allowing ideal spin injection and detection. 

The Landauer-Buttiker formalism is used for the numerical conductance calculation and 

is used to evaluate matrix elements of related Greens' functions. 
The following parameters are used; m* = 0.04 x m electron , where m electron is the free electron 
mass, and the Fermi energy Ep = 0.103 eV, which, in the absence of a and V c , amounts 
to the Fermi wavelength 2Tx/kp = 191 A and the electron density of n s = 1.72 x 10 12 cm -2 
in the 2DEG. For the TB approximation, we choose a = (27r//cj?)/9.95. Thus the hopping 
energy t = h 2 /2m*a 2 becomes 0.259 eV. 

Below the conductance is first evaluated in two windows of a, (0~10)xl0 -12 eVm and 



the method in Ref. 



(45~55) x 10 12 eVm. In various candidate systems 



InAs/GaSb [11], InGaAs/InP lj, and InAs/AlSb 



or the SFET such as InGaAs/InAlAs 



a, 



is| , a typically falls within the first win- 
dow. A recent calculation [23J also reported that for III-V quantum wires, a is typically 
in this weak coupling regime. Search for systems with higher a is in progress and ef- 



fects of strong RSO coupling on the energy spectrum are addressed theoretically |24j. For 
CdTe/HgTe/CdTe, for instance, values up to ~ 45 x 10~ 12 eVm has been reported jl^ . 
which motivates the second window. For a SFET with relatively short L, the conductance 
is evaluated in a much wider window of a, (0~100)xl0~ 12 eVm. 

Figure 121 shows the evolution of the conductance G as a function of a for the SFET 
with 3, 5, 10, and 20 channels, respectively. For the range of a examined below, the 
number of channels is related to the width w of the 2DEG via the relation j^j N c ^ = 
lnt(kpw/7[), where Int(x) denotes the largest integer not exceeding x. The length L of 
the 2DEG is assumed to be 1.44 /xm for (a) and (b), and 0.8 /im (w 0.556 x 1.44 /jm) for 
(c), and 0.1 x 1.44 /im for (d), respectively. In Fig. EJa), the SFET is operated in the a 
range from to 10 x 10~ 12 eVm. To estimate the spin precession angle, we use the formula 
2m*aL/h 2 . Though this formula is derived in the single-channel limit, it may still be useful 



to estimate the spin precession in multi-channel SFETs. Then the average spin precession 
angle, estimated from the formula 2m*aL/h 2 with the average a = 5.0 x 10~ 12 eVm, is 
around 1.2 x 2ir, and the variation of the spin precession angle within the specified range 
is ±1.2 x 2tt. Note that the SFET exhibits an almost ideal conductance modulation in the 
sense that the maximum of G reaches almost {e 2 /h)N c ^ and the minimum of G reaches 
almost zero. The conductance modulation behavior becomes, however, less ideal with the 
increase of iV c h. For = 20, for instance, Fig.EJa) already shows some deviation from the 
ideal behavior. Figure Efa) shows, as a function of N ch , the ratio between G max and G m i n , 
where G max and G m ; n are the maximum and minimum values of G in the interval 1.8 x 10 -12 
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FIG. 2: (Color online) Conductance G of a ballistic multi-channel SFET as a function of the RSO 
coupling parameter a at the temperature T=0. Ideal spin injection and detection are assumed. 
The following parameters are used; the effective mass m* = 0.04 x m e i ectr on, where m e i e ctron is the 
free electron mass, and the electron density n s = 1.72 x 10 12 cm~ 2 . The length L of the 2DEG is 
1.44 ^m for (a) and (b), 0.8 /jm (« 0.556 x 1.44 /im) for (c), and 0.1 x 1.44 /jm for (d). 



7 



CD 

>< 
to 



i 

E 
CD 



CD 

x 
to 

E 

CD 



10° 
10 4 
10 2 

1 

10 6 

1(f 
10 2 

1 

10 6 
10 4 
10 2 




H 1 1 1 h 



20 



40 



60 



(b) 







(c) 



80 100 



N, 



'ch 



FIG. 3: (Color online) The ratio between the maximum conductance G max and the minimum 
conductance G m j n as a function of the number of channels The panels (a), (b), and (c) 

correspond to the situations in Fig. Ufa), (b), and (d), respectively. The parameters are the same 
as those for Fig. |2J 

eVm < a < 4.5 x 10~ 12 eVm in Fig. 12(a). It clearly shows the decay of the ratio with 

In Fig. EJb) (again with L = 1.44 /xm), the SFET is operated in the a range (45 ~ 
55) x 10 -12 eVm. Note that the width of the a range is the same as that in Fig.^Ja). Since 
a is larger than the values in Fig. E^a), the spin precession angle will be larger. The average 
spin precession angle in this range, estimated from the formula 2m*aL/h 2 with the average 
a = 50 x 10~ 12 eVm, is around 12 x 2tt, about 10 times larger than in Fig. 12(a). On the 
other hand, the variation of the spin precession angle within the specified range is about 
±1.2 x 27r, same as in Fig. |2fa). Note that the conductance modulation behavior is now 
much less ideal; the maximum and minimum of G deviate considerably from (e 2 //i)iV c h and 
0, respectively. Moreover the conductance oscillation with a becomes somewhat irregular. 
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Figures b) shows, as a function of N ca , the ratio between G max and G m i n for the SFET 
in Fig. I2fb). Here G max and G m i n are evaluated in the a range (48.5 ~ 51.5) x 1CT 12 eVm. 
Note that the decay of the ratio with increasing N ca is much faster than that in Fig. [5^ a). 
For small N cri around 5, the ratio G max /G m i n already drops to of order unity. 

One intuitive explanation for the much less ideal behavior in Figs. I2fb) and|Hfb), com- 
pared to Figs. |2fa) and Of a), is fluctuations of the spin precession rates from channel to 
channel. For example, if the spin precession angles per length deviate from 2m* a /h 2 and 
the deviations fluctuate from channel to channel, the conductance modulation may become 
less ideal. To test this possibility, we examine a SFET with shorter L since effects of the 
channel-to-channel fluctuation are expected to be weakened with decreasing L. In Fig. EJc), 
the length L of the 2DEG is reduced to 0.8 /im (» 0.556 x 1.44 /mi). When the SFET is 
operated in the same a range as in Fig. I2fb), the average spin precession is estimated to 
6.7 x 27r, which is about 44% smaller than that in Fig. Efb), and the variation of the spin 
precession angle within the specified range is estimated to ±0.7 x 2n. Figure|2^c) shows that 
the conductance modulation is still far from ideal despite the 44% reduction in the average 
spin precession angle. 

In Fig. I2fd), the length L of the 2DEG is further reduced to 0.1 x 1.44 pan., one tenth of 
that in Fig.|2^a). When the SFET is operated in the range < a < 100 x 10~ 12 eVm, which 
is 10 times larger than in Fig. [2Ja>), the average spin precession angle and the variation of 
the angle within the specified a range are exactly same as those in Fig. CD^a). Thus the 
situations in Figs. Efa) and (d) are identical as far as the expression 2m*aL/h 2 for the 
spin precession angle is concerned. The numerical calculation result in Fig. Efd) indicates 
however that while the conductance modulation behavior is somewhat improved compared 
to those in Figs. [2Kb) and (c), it is still less ideal than that in Fig. Ufa). Figure Efc) shows 
the iV ch - dependence of the ratio G max /G min for the situation in Fig. Efd), where G max and 
G m m are evaluated in the interval 18 x 10~ 12 eV < a < 45 x 10~ 12 eVm. Note that the decay 
of the ratio is still much faster than in Fig. Efa). This result leads one to conclude that the 
channel-to-channel fluctuations of the spin precession angles are not the main reason for the 
deviations from the ideal behavior. 
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B. Transport channels 

In order to gain an insight into the numerical results in Figs. El and H3 it is useful to 
examine energy eigenfunctions and eigenenergies of H 2 r> [Eq. Let ^f(x, y) = e lkxX ifj(y) 
denote a eigen wavefunction with energy E. Due to the translational symmetry along the 
x-direction within the 2DEG, it can be expressed as a superposition of four plane waves 
^ = X/j=i c i^^v where the plane wave ^^{x^y) = e lkj ' r Uj(\tj) is an eigenstate of H 2 d in 
the absence of V c (y). Here all four wavevectors k, = (k x j, k y j) share the same longitudinal 
component k x j = k x and their transverse components k y j are determined from the relation 
h 2 {k 2 - + k 2 -)/2m* + (—iy<y(kl j + k 2 -) 1 / 2 = E, which is the energy-dispersion relation in 
the absence of V c (y) (see Fig. EJ). The spinor Uj(k) is given by 



where <^(k) is defined by cosv?(k) = k x j{k 2 x + k 2 ) 1 / 2 and sin^(k) = k y /{k 2 + k 2 ) 1 / 2 . Then 
the four constraints from the boundary conditions ty(x,0) = ty(x,w) = fix the four 
coefficients Cj, and the exact energy-dispersion relations and the eigen wavefunctions can 
be obtained. This procedure can be simplified further by using additional symmetries of 
if 2D- See Appendix IB1 for further details on the symmetries, and Appendix IU1 for the use 
of the symmetries for the evaluation of the exact energy-dispersion relation and the eigen 
wavefunctions. 

The red dotted lines in Fig. El denote the exact energy- dispersion relation for various 
situations. In comparison, the perturbatively obtained energy-dispersion relation (up to the 



FIG. 4: (Color online) Four solutions k y>j (j = 1,2,3,4) of h 2 {k 2 x +k 2 yj ) /2m* + {-l)ia{k 2 x +k 2 yj ) 1 1 2 = 
E, which is the energy dispersion relation in the absence of V c (y). Note that for certain k x (denoted 
by the red dashed line), two solutions become purely imaginary numbers. 
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second order in a) , 



E 



n,i 



h 2 k 2 

x +£ n +(-lYak x - 



m*a 2 



2m * ■ -•- ■ v -/ 2 h 2 ' ^ 

is also depicted (blue solid lines) in Fig. El Here £ n = h 2 k 2 y /2m* is the energy of the n-th 
(n— 0,1,2,- ■ •) excitation mode in the transverse direction, k n , y = (n + l)ir/w, and i — 1 or 
2. Equation (JBJ) is obtained by treating the RSO coupling term in H 2 r> as a perturbation. 
Further details on the calculation of the energy dispersion relation via the second-order 
degenerate perturbation theory are given in Appendix [Dl 

Figure El shows the local spin directions of the selected eigen wavefunctions at the Fermi 



0.398 




M 1 ° m ) 



/c x (10 8 m" 1 ) 



FIG. 5: (Color online) The energy dispersion relations in a Rashba spin-orbit (RSO) coupled system 
with hard wall confinement. Red dotted lines represent the exact dispersion relations, while blue 
solid lines the perturbative results, Eq. jSJ). The two types of the lines overlap almost perfectly in 
(a) and (b). The RSO coupling parameter a and the number of channels are a = 8.3 x 10~ 12 
eVm, N ch = 3 in (a), a = 8.3 x 1(T 12 eVm, N ch = 10 in (b), a = 49.8 x KT 12 eVm, iV ch = 3 in 
(c), and a = 49.8 x 1CP 12 eVm, iV c h = 10 in (d), respectively. The width w is 32.7 nm for (a) and 
(c), and 100.0 nm for (b) and (d). The energy is given in units of t = 0.259 eV. The horizontal 
dashed lines denote the Fermi energy Ep = 0.398i = 0.103 eV. 
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level Ep. The width w and the RSO parameter a for Figs. Efa), (b), (c), and (d) are the 
same as those for Figs. 0a), (b), (c), and (d), respectively. The local spin direction of an 
energy eigen wavefunction is independent of x due to the translational symmetry of H 2 d 
along the x-direction. Also it can be shown that it is always perpendicular to the x-axis 
and remains within the yz-plane since H 2 b remains invariant under the symmetry operator 
QU x V spm ('k, n) (see Appendix El for further details), where G is the time-reversal operator, 
U x is the mirror reflection operator with respect to the yz-plsnae in the orbital space, and 
X> spin (x, 7r) is the rotation operator with respect to the x-axis by the angle ir in the spin 
space. In each upper panel in Figs. Efa) , (b), (c), and (d), the arrows indicate the local 
spin direction as a function of y for six selected exact eigen wavef unctions. The length of 
the arrows is proportional to ^^(x, y)^/(x, y), which is independent of x. The six states 
in Figs. Efa) and (c) are the six eigenstates in Figs. 02a) and (c) at Ep = 0.398t = 0.103 
eV with increasing order of k x > . On the other hand, the six states in Figs. Efb) and 
(d) are the eigenstates in Figs. Efb) and (d) at Ep with the smallest, 2nd smallest, 7th 
smallest, 8th smallest, 13th smallest, and 14th smallest k x s (ordered according to the exact 
k x values) . The lower panels in Figs. Efa), (b), (c), and (d) show the local spin angle 9 as 
a function of y, where 9 is —tt/2, 0, tt/2, ±tt when the local spin direction is pointing — y, 
+z, +y, —z axis, respectively. Here 9 is related to the local spin direction via the spinor 
representation (cos |xt + i sin ^Xi) °f a local spin direction that lies within the yz-plane, and 
Xi — y I j , Xl = \ are ^ ne s Pi nors pointing along the (±z)-directions. In each lower panel 
of Figs. EJa), (b), (c), and (d), there are three subpanels, each of which shows the profiles 
of 9(y) , s for two selected energy eigen wavef unctions; #'s in the left /middle/right subpanel 
correspond to the two leftmost /central/rightmost spin textures in the upper panel. In the 
lower panels, blue dotted lines represent the exact results (see AppendixOfor further details) 
and red solid lines are from the perturbatively obtained eigen wavefunctions calculated up 
to the second order in a (see Appendix O for further details). 

In Figs. El and El two representative values of a, 8.3 x 10~ 12 eVm [(a) and (b)] and 
49.8 x 10~ 12 eVm [(c) and (d)], are used. Figures Ef a) and Efa) show the result for the 
smaller a value and iV c h = 3. Note that the perturbative result for the energy dispersion 
relation is in excellent agreement with the exact result. Equation (JHJ) from the perturbation 
theory indicates that the energy dispersion relations for all channel n are identical to each 
other up to the parallel shift by S n and remain parabolic. Moveover Fig. Ef a) indicates that 
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FIG. 6: (Color online) Spin texture along the transverse direction for selected states at Ep. The 
parameters m*, a, w and in (a), (b), (c) and (d) are the same as those in Figs. Eta), (b), (c) 
and (d), respectively. In the lower panels, the blue dotted lines represent the exact results while 
the red solid lines the perturbative results. Inset in the upper panel of (c) shows the magnified 
structure near the points where y/^(x, yf¥{x, y) almost vanishes. 
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the spin direction is essentially equal to the (+y) or (— y)-direction, same as in the ideal 
single-channel SFET. Thus the SFET for the situation in Fig. 02a) andlf^a) should show an 
almost ideal behavior, which is indeed the demonstrated in Figs. |5(a) for N ch = 3. 

The channel-to-channel fluctuations of the spin precession angle should be negligible. 

Figures EJb) and EJb) show the result for a = 8.3 x 10 -12 eVm and = 10. The 
exact energy dispersion relations are again excellently fitted by the perturbative results 
[Eq. (JHJ)], which are purely parabolic. The local spin directions, on the other hand, show 
noticeable deviations from the ideal (±y)-directions. The lower subpanels in Fig. EJb) show 
that near the "nodal" points, where the length of the spin arrows in the upper panel of 
Fig. Efb) almost vanishes, 9(y) changes rapidly by 2tt, forming "phase-slip" -like structures. 
The phase- slip-like structure near the nodal points is more clearly visible in the magnified 
plot in the inset of Fig. Efc). Superimposed over the phase-slip-like structures, there is also 
a gradual drift of 6(y) with the increase of y. To understand the origin of the deviations, the 
eigen wavefunctions are perturbatively evaluated (see Appendix |Dj). Up to the first order 
in a, the perturbative calculation for the reduced wavefunction if>(y) [&(x, y) = e lkxX ip(y)} 
results in 

where n and i are defined in the same way as in Eq. (jHJ), and i = 2 for i — 1 and i = 1 
for 2 = 2. Here the zeroth order reduced wavefunctions are ip^iiv) = 4> n (y) • (xt + *Xi)/V2 
and ipn\{v) = 4>n{y) ■ (xt —iXl)/V2 respectively, where (/) n (y) = y/2/w sin k n , y y- From these 
expressions, one obtains simple expressions for 6, 

7T 9?77*rv 

e nti =i(y) = 2 - -jr^ - w ^ + o(a 2 ), (io) 

9 n , t=2 (y) = -\- 2 -^{y-w/2) + 0{a 2 ), 



which has been reported previously 



25| . Note that while the zeroth order terms in a are 



^-independent constants, ±vr/2, corresponding to the ideal (±y)-directions, the first order 
terms in a are proportional to y — w/2, explaining the gradual drift of 6 with y. In order 
to gain an insight into the phase-slip-like structure, we next examine the second order 
corrections to ip(y). Since the first order correction in Eq. Q identically vanishes at the 
nodal points (sin k n ^y = 0) of ^°], the second order corrections are the first nonvanishing 
terms near those nodal points. Moreover since if) = is a topologically singular point for 
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the local spin angle 6, even a small and smooth corrections to if) can result in a rapid change 
of 9 by 27T. The local spin angles evaluated from the second order perturbation theory are 
shown (red solid lines) in the lower subpanels in Fig.EJb), indeed reproducing the phase-slip- 
like structure. The perturbative results are in reasonable agreement with the exact results. 
The expressions for the second order corrections to ip are rather lengthy [Eqs. (|D14|) and 
(|D15|) ] and given in the Appendix [D] only. We remark that the gradual drift and the phase- 
slip-like structures are also present in Fig. EJa), though they are much weaker. Since they 
can be explained within the perturbation theory, the strength of their effects is expected 
to depend on the magnitude of the small expansion parameter of the perturbation theory, 
which is 2m*aw/h 2 according to Eqs. Q, ()D14|) . and (|D15|) . Then recalling that the product 
2m*aw/h 2 is about 0.28 for the 3 channel system and about 0.86 for the 10 channel system, 
the difference between Figs. Efa) and (b) can be explained. This analysis also implies that 
the deviations of the local spin configuration from the ideal spin directions will be magnified 
with the increase of (oc w). Since these deviation will certainly cause the conductance 
to deviate from the ideal SFET behavior based on the ideal spin directions, this analysis 
also provides an explanation for the initial decay of the conductance modulation ratio with 
the increase of N ch demonstrated in Figs. Ufa) and Ota). 

Figures Efc) and^c) show the result for the larger a = 49.8 x 10~ 12 eVm and N c ^ = 3 
with the product 2m*aw/h 2 ~ 1.7, which is about two times larger than in Figs. Efb) 
andlH^b). Figure Efc) shows that now some energy dispersion relations noticeably deviate 
from the parabolic behaviors. While many subbands are still well fitted by the parabolic 
perturbative results, the two subbands, which are predicted by the perturbation calculation 
to cross near k x « 2.5 x 10 8 m _1 , are affected by the subband mixing and the avoided 
crossing structure appears, which goes beyond the perturbation calculation in Appendix [DJ 
At E F , the states with the second and third largest k x 's are affected by the avoided crossing 
and their dispersions deviate from the parabolic dependence, modifying the spin precession 
rates of the involved channels as demonstrated recently |26j. The avoided crossing also 
affects the local spin angle. For those two states affected by the avoided crossing, the exact 
local spin angle profiles deviate considerably from the perturbative results [lower subpanels 
of Fig. Efc)]. For the other four states, the exact local spin profiles are reasonably well 
fitted by the perturbative results. Note that the gradual drift of 9 and the phase-slip-like 
structure are more evident than in Fig. Efb), which is natural in view of the about two- 
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fold increase of the product 2m*aw/h 2 . Based on this analysis, we attribute the reduced 
conductance modulation ratio in Figs. I2jb) and Efb) (even for small A" c h such as 3) jointly 
to the perturbative spin configuration change and the avoided crossing. We remark that 
the avoided crossing in fact occurs in Fig. E(b) as well, though not clearly visible in the 
figure since its effects on the energy dispersions are restricted to rather narrow ranges of k x . 
Furthermore the avoided crossing appears far below Ep, so that its effects on the conductance 
is negligible. 

Figures Efd) and Efd) show the result for a = 49.8 x 10~ 12 eVm and = 10 with 
2m*aw/h 2 ~ 5.2. It turns out that most subbands are affected by the avoided crossing 
and the agreement between the perturbative and exact dispersion relations is poor. The 
agreement between the perturbative and exact spin configuration is again poor for many 
states at Ep. The abundance of the avoided crossing explains the strongly suppressed 
conductance modulation ratio of the SFET operated with larger a [Figs. Efb) and (c)]. 



III. COHERENT SFET 

The electronic coherence may give rise to interesting effects in mesoscopic systems j^oi I27 1 



For a single channel system, it has been predicted [15J, [16|, |17|, [18| that due to the electron 
coherence, the conductance behavior of a SFET deviates from the conventional sinusoidal 
dependence cos 2 (m*aL/h 2 ). In this section, we investigate effects of the coherence on multi- 
channel SFETs. In particular, we focus on situations with N c ^ < 10 and a < 10 x 10~ 12 
eVm since for larger and a, the coherence effects will be mixed with other effects such as 
nonideal spin configuration and avoided crossing, making clear identification of the coherence 
effects difficult. 

Hints of the coherence effects appear already at Fig. |2fa). A close look reveals the 
deviation of the line profile G(a) from the sinusoidal behavior oc cos 2 (m* a L/h 2 ). The 
evolution of G(a) with the variation of the Fermi energy Ep illustrates the deviation more 
clearly. Figure [7)^a) shows G[a) of a SFET with A" c h = 3 for a number of Ep's. 5Ep's in 
Fig. Ufa) are measured with respect to the reference energy 0.103 eV=0.398t. While G(a) 
is not sensitive to Ep near the conductance maximum (near a = 4.2 x 10~ 12 eVm) and the 
minimum (near a = 2.2 x 10~ 12 eVm and 6.2 x 10~ 12 eVm), the profiles of G(a) including the 
full-width-at-half-maximum (FWHM) exhibit considerable dependence on Ep. Moreover 
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FIG. 7: (Color online) (a) Variations of the conductance profile with Ep for a SFET with = 3 
and L = 1.44 /im at the temperature T = 0. 5Ep J s are measured with respect to the reference 
energy 0.398i = 0.103 eV. Inset : The finite temperature conductance G(a) at kgT = IOHvf/L and 
bEpjt = 0.209 where vp = H(Kq + K\ + K2)/3m*. Note that the numerical result is in excellent 
agreement with the analytic high temperature expression 28]. (b) The red solid line represents the 
numerically calculated G(a* = 3.94 x 10 -12 eVm) as a function of Ep at T = 0. Here a* [vertical 
dotted line in (a)] is close to the conductance maximum point. Thus dips in the G(a*) graph 
imply the appearance of the nested peak structure at the corresponding Ep. For comparison, the 
prediction of the formula [Eq. (|12j) ] is also shown (blue dashed line), (c) sin 2 (K n L) for n =0, 1, 
and 2. The comparison between (b) and (c) indicates that the dip positions of G(a*) coincide with 
the points where one of sin 2 (i ; C n L)'s becomes zero. 

near certain special Ep's such as 5Ep/t = 0.209, G(a) shows a nested peak structure, 
where a small peak appears on top of a bigger background peak. These results are in 
clear contrast to the result in Ref. Q], which reports the absence of the ^-dependence 
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when 2m*aw/h 2 <C 2tt. While our conductance calculation takes into account all three 
parts (2DEG, injector, and collector) quantum mechanically, only the 2DEG is quantum 
mechanically treated in Ref. We believe that this is responsible for the difference. 

In order to examine the connection between the electronic coherence and the nested peak 
structure, we first recall that in a single channel SFET, the conductance calculation jl8j |. 
which takes into account the electronic coherence, results in 



where 4> a = 2m*aL/h 2 , K = [{k x ) i= i + (k x )i=2]/2, and {k x )i=\ and {k x ) i= 2 (> 0) are the two 
longitudinal wavevectors at Ep. The deviation of Eq. from the sinusoidal behavior arises 
from the fact that the injector and collector behave as reflecting walls for the electron spin 
component antiparallel to the favored spin direction in the injector and collector. Thus the 
unfavored spin component may be reflected many times by the ideal injector and collector 
until it acquires, via the spin precession by the RSO coupling, some component parallel 
to the favored spin component and is transmitted to the injector or collector. Thus in 
addition to the direct contribution that does not undergo any reflection, there are n-th 
order contributions to the electron transmission that undergo the reflection 2n times until 
being transmitted to the collector. In this sense, the coherent SFET is analogous (though 
not identical) to the Fabry-Perot interferometer. It is straightforward to verify that the 
amplitude of the (n + l)-th contribution is smaller than that of the n-th contribution by 
the factor — exp(i2KL) sin 2 (0 Q /2) for n = and by the factor exp(i2KL) cos 2 (</> Q /2) for 
n > 1. Here the zeroth (n = 0) contribution denotes the direct contribution. Note that 
the interference between the direct contribution and the multiply reflected contributions is 
destructive (constructive) when the phase exp(i2KL) is +1 (—1). This explains the factor 
sin 2 (KL) in Eq. (fTT|) . which predicts the suppression (enhancement) of the FWHM of the 
profile G(a) when sin 2 (KL) « (1). Since K depends on Ep, this K dependence also 
implies the Ep dependence. To examine the coherence effects in a multi-channel SFET, we 
use an intuitive generalization of Eq. to a multi-channel system; 



where K n = [(k x ) ni i=i+(h x ) ni i=2\/2, and {k x ) n ^ (> 0) is the longitudinal wavevector at Ep for 
the subband (n, i) (see Fig. EJ). Equation (fT2^) is expected to be a reasonable approximation 



G 



4 cos 2 \ sin 2 (KL) 



(11) 



h sin 4 if + 4 cos 2 ^ sin 2 (XL) 




(12) 
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when the inter-channel mixing is weak and each channel is close to an ideal ID channel, 
which are indeed the case for a < 10 x 10~ 12 eVm and iV c h < 10 as demonstrated in 
Sec. H2 In Fig- ITJb), the prediction of Eq. (fP2j) (blue dashed line) for G(a*) as a function 
of Ep is compared with the numerical TB conductance result (red solid line), where a* = 
3.94 x 10~ 12 eVm [vertical dashed line in Fig. [7(a)]. Since a* is close to the position where the 
conductance retains its maximal value (e 2 /h)N c \ 1 , the dips in G(a*) imply the formation of 
the nested peak structure. To compare with the TB calculation result, K n and L in Eq. (|T2"j) 
are expressed in terms of the TB parameters. To evaluate K n = [{k x ) n ^ i= i + {k x ) n ^ = 2\/2, we 
use the relation Ep = 2t {2 — cos[(k x ) nti a] — cos[(n + l)7ra/u>]} + (— l)*(a/a) sm[(k x ) nti a} — 
m*a 2 /2h 2 , which amounts to the TB approximation of Eq. (jSJ). Here w = a(N y + 1). To 
evaluate L, we use L = a(N x + 1 + (3), where (3 = 0.6 is an ad hoc fitting parameter much 
smaller than N x = 749. Note that Eq. ()12J) produces dips (blue dashed line) and moreover 
the dip positions predicted by Eq. (|T2"j) are in good agreement with the TB conductance 
calculation result (red solid line). Thus Eq. fT2|) reproduces the nested peak structure. We 
remark that when f3 = is used, the dips in the blue dashed line shift to the right by 
0.0001 ~ 0.0003/f: only (not shown), not affecting the agreement significantly. 

The agreement provides a simple insight into the origin of the nested peak structure. A 
key observation is that while all g^s in Eq. (I12j) produce the conductance peaks at the same 
a, the peak widths may fluctuate considerably from channel to channel. The channel-to- 
channel peak width fluctuation arises from the dependence of the peak width on K n . In 
particular when sin 2 (K n L) w 0, the peak produced by g n becomes very narrow. Thus when 
sin 2 (i^ ra L) is close to zero for one particular n, one very narrow peak for that particular n is 
superposed with the other broad peaks, generating the nested peak structure. Figure [7(c) 
plots sin 2 (fT n L)'s as a function of Ep. Note that whenever one of sin 2 (i^ n L) vanishes, a dip 
appears in Fig.UJb). We remark that though the nested peak structure is illustrated mainly 
for iV c h = 3, it appears for higher iV c h as well. An example for A^h = 20 is visible in Fig. [2(a) 
near a = 4.3 x 10~ 12 eVm. For sufficiently large iV c h, however, the nested peak structure is 
expected to vanish. The nested peak structure was not found in a two-dimensional coherent 
SFET j^l, which amounts to the N c ^ —>■ oo limit. 

The electronic coherence can give rise to another interesting effect when an external 
magnetic field B is applied, providing an additional source of the spin precession. Here 
we consider B = B\\y applied along the j/-axis, so that the spin precession axis by B is 
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FIG. 8: (Color online) Conductance G as a function of the RSO coupling parameter a for (a) 
iVch = 3 and (b) = 10. The solid lines in (a) and (b) show the analytic results obtained using 
Eq. (JT3J) for a specific field strength 2g*fj, B B\\/t = 1.33 x 1CT 3 and 1.38 x 1CT 3 , respectively. Note 
that the analytic results are in good agreement with the numerical results. 



parallel to that due to the RSO coupling in the ideal ID situation. The spin precessions 
by B» and the RSO coupling differ in one crucial respect; When the electron reverses its 
motion, the direction of the spin precession by the RSO coupling is reversed while that 
of the spin precession by B\\ is not. For a single channel SFET, it is illustrated 3] that 
this difference, combined with the coherent Fabry-Perot-like interference, gives rise to the 
splitting of the conductance peaks upon the application of Bn . Figure |H1 shows that the 
magnetic-field-induced peak splitting occurs also for multi-channel SFETs. The numerical 
TB conductance calculation results are compared with the expression, 

2 iVch-1 



2 N ch -1 



n=0 



e 

77 



^ 4 (cos 2 4f - sin 2 sin 2 K n L + 4 sin 2 4f sin 2 *f 



n=Q 
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2 4f + sin 2 %^) 2 + 4(cos 2 4f - sin 2 sin 2 K n L 



which is an intuitive extension of the single channel formula in Ref. 



(13) 



18| . Here (f)B, n denotes 
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the spin precession angle due to B\\. For weak fields, it is given by 2g* /jlbB\\/ (hvF, n / L) , where 
g* is the effective Lande-g factor, is the Bohr magneton, and vf,u = hK n /m* is the group 
velocity of electrons with energy Ep in the n-th channel Note that the predictions of 
Eq. (J13|) are in reasonable agreement with the numerical TB results both for N c ^ = 3 and 
10. 



IV. DISCUSSION 



An ideal spin injector /collector has been assumed in the preceding sections in order to fo- 
cus on the multi-channel effects. For a certain material to be an ideal spin injector/collector, 
(i) its spin polarization should be 100% and (ii) its effective mass and Fermi wavelength 
should match those for the 2DEG. While the condition (i) is satisfied in the so-called half- 
metals (jyj and the condition (ii) may be satisfied in ferromagnetic semiconductors or diluted 
magnetic semiconductors j^, we are not aware of any material that satisfies the both con- 
ditions. In Sec. IIV A| we thus address briefly effects of a nonideal spin injector /collector on 
the multi-channel SFET. 



A. Nonideal spin injector/collector 

One of most explored and representative choices for a spin injector /collector is conven- 
tional metallic ferromagnets, whose spin polarization can be as high as ~ 50%. Use of 
conventional metallic ferromagnets requires a special care. Recent theories proposed 
that the spin injecton/detection rate can be greatly improved by introducing a thin insu- 
lator between a conventional ferromagnet and the 2DEG while without the insulator, the 
spin injection/detection rate was demonstrated 5] to be below 1% due to the so-called 
conductance mismatch problem The physics behind the insulator "prescription" for 
the conductance mismatch problem is that the tunneling through the insulating barrier be- 
comes spin-dependent Q in the fe™magne W nsulato r -2DEG geometry even though the 
insulator is nonmagnetic. Subsequent experiments |24J using oxide insulating barriers have 
reported the enhanced spin polarization of 2-30% at room temperatures. Schottky tunnel 
barriers were also demonstrated to be effective and the spin polarization of 2-30% has been 



reported |35J] at low temperatures. 
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Here we calculate numerically the conductance of the multi-channel SFET equipped with 
the nonideal spin injector /collector that consists of a conventional metallic ferromagnet-thin 
insulator-2DEG hybrid structure. One possible way to model the spin- dependent tunnel- 
ing in the ferromagnet-insulator-2DEG structure is to explicitly take into account all three 
materials, ferromagnet, insulator, and 2DEG (Appendix [E]) . This approach is however 
too demanding for numerical calculations since the Fermi wavelength in the ferromagnet is 
about two orders of magnitude smaller than that in the 2DEG. We thus adopt a simplified 
phenomenological description, where all effects of the ferromagnet (including not only its 
spin-dependent density of states but also the effective mass and the Fermi wavelength differ- 
ences from those of the 2DEG) are taken into account by phenomenological spin- dependent 
hopping parameters for the tunneling through the insulator. To be specific, the TB Hamil- 
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FIG. 9: (Color online) The conductance G of a ballistic SFET with the nonideal injector /collector 
described by the modified #™ pling [Eq. dJ)] and ff?f /col . iVch = 3 in (a) and N ch = 10 in (b). The 
length L of the 2DEG is 1.44 in both (a) and (b). Other parameters are the same as those in 
Fig. [3 The spin polarization of the injected current for the ferromagnet-insulator-2DEG structure 
is assumed to be about 30%. 
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FIG. 10: (Color online) The conductance G as a function of the Fermi energy Ep with fixed a, 
a m in = 2.07 x 10~ 12 eVm and a max = 4.13 x 10~ 12 eVm, for N ch = 3 in (a) and N ch = 10 in (b). 
The peak positions of G for a max are marked by arrows. Note that there are three closely spaced 
peaks near E F = 0.397t. The length L of the 2DEG is 1.44 /xm. 



tonian [Eq. 
as follows: 



for the coupling between the 2DEG and the injector /collector is modified 



rrTB 

coupling 



" Z) X) *«« (''I-/- - + C N*+l,j,s x C N x j,s* + H.C.J , (14) 

where t Sl is a spin-dependent hopping parameter that takes into account not only the spin- 
dependent tunneling through the insulator but also other effects from the effective mass and 
the Fermi wavelength differences between the ferromagnet and the 2DEG. Since all effects of 
the ferromagnet are already taken into account in t Sx [Eq. (|14j) ]. the injector /collector may 
be modelled as a nonmagnetic semiconducting lead, which can be described by modifying 
if™ col [Eqs. (HJ) and (j3J)] so that similar contributions for s x = —1 are also included. Al- 
though this description adopts a considerably simplified picture of the real system, it is still 
expected to capture the main physics of the spin injection/detection since, as demonstrated 
m Refs. Q and Appendix |Ej the spin- dependent tunneling through the insulator is the 
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FIG. 11: (Color online) (a) The conductance G as a function of Ep for the SFET with two nonideal 
spin injector/collector that contains a thin insulating barrier. Here = 10 and a = 8.3 x 10~ 12 
eVm. For the curves labelled as the strong and weak barriers, the hopping parameters [Eq. (|14|) ] 
are (t Sx =i = O.li, t Sx= -i = 0.073t) and (t Sx= i = 3 x O.lt, t Sx= -i = 3 x 0.073i), respectively. Inset in 
(a) : G{Ep) for = 3. Note that the positions of the large resonance peaks nearly coincide with 
those of the nested peaks in Fig. Efb). (b) The conductance G as a function of a for the strong 
and weak barriers when = 10 and Ep = 0.416i. The length L of the 2DEG is 1.44 /xm. 

main source of the sizable spin injection/detection rate (~ 30%) while the bulk properties of 
the ferromagnet generate much weaker contributions only. We choose the value of t Sx to be 
O.lt for s x = 1 and 0.073t for s x = —1, respectively. The resulting spin polarization of the 
injected current estimated by (tf =1 — tl _ _i )/{tl _ 1 + if = _x) is then 30%, which is similar 
to the values reported in experiments [34|. 

For iVch = 3 and Ep = 0.398i [solid line in Fig. EH^a)], G clearly shows the modulation 
by a though the value of G is considerably reduced below 0.001 x (e 2 /h) due to the insu- 
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FIG. 12: (Color online) The conductance G of a ballistic SFET with the nonideal injector /collector 
as a function of the RSO coupling parameter a at T = 0. The length L of the 2DEG is 0.1 x 1.44 /iin 
and Ep = 0.398i = 0.103 eV. 

lating barriers and the ratio G max / G m i n ~ 1.7 is also suppressed due to the reduced spin 
injection/detection rate of 30%. Thus for this particular situation, G shows rather expected 
behaviors. For other situations, somewhat unexpected behaviors also appear. For N c ^ = 3 
and E F = 0.416t [dashed line in Fig. Eta)], the "periodic" modulation disappears and in- 
stead a large peak appears at a ~ 2.1 x 10~ 12 eVm with the maximum conductance of 
~ 1.8 x (e 2 /h) (the peak height is scaled down by the factor of 5000 to fit in the vertical 
scale of the graph). For iV c h = 10 and Ep = 0.416t [dashed line in Fig. E^b)], a large peak 
(near a ~ 8.9 x 10~ 12 eVm) coexists with the G modulation in the lower a range. For 
iV c h = 10 and Ep = 0.398t [solid line in Fig. Efb)], the G modulation is superimposed on a 
background slope. 

To investigate how frequently such peaks appear, G is evaluated as a function of E F for 
two representative values of a, a min = 2.07 x 10~ 12 eVm and a max = 4.13 x 10 -12 eVm, 
which correspond to the positions of the local minimum and maximum of G at Ep = 0.398t 
for iV c h = 3. For A^h = 3 [Fig. ITuTa)]. two peaks appear in the depicted Ep range and 
for iV c h = 10 [Fig. ITW b)]. seven peaks appear. Thus the occurrence of the peaks becomes 
more frequent with the increase of N c ^. The calculation for iV c h=3, 10, and 20 over a 
much wider range of Ep (0.360t < Ep < 0.440t, not shown) verifies this trend. Note 
that the peak positions for a m ; n and a max are somewhat different, implying that the peak 
positions shift with a. Considering that the peak heights (~ e 2 /h) are three or four orders of 
magnitudes larger than the modulation amplitudes in Fig. |Hl even tails of the peaks can be 
large "perturbations" to the modulation and thus the a-dependence of the peak positions, 
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FIG. 13: (Color online) The conductance G for several impurity-configurations in the 2DEG as 
a function of the RSO coupling parameter a for iV c h = 3 in (a) and iV c h = 10 in (b). The seed 
numbers represent different impurity configurations, all of which have the same mean free path 
I = L/3. The other parameters are the same as in Fig. Efa) and (b), respectively. The length L of 
the 2DEG is 1.44 /xm. 

though weak, can affect the modulation considerably. As remarked above, this perturbation 
appears more frequently for larger 

To find the origin of the large peaks, their dependence on the strengths of the hopping pa- 
rameters t Sx is examined. As demonstrated in Fig. ^2 t ne peak widths for the large hopping 
parameters {t Sx= \ = 3 x O.lt, i Sa .=_i = 3 x 0.073t) are larger than those for the small hopping 
parameters (t Sx= i = 0.1t, t Sa .=-i = 0.073£). This difference in the peak widths implies that the 
large peaks are resonances due to the electron confinement effects of the insulators. Then 
at the energies where the resonances appear, the phase acquired by electrons during their 
motion from the injector to the collector should be an integer multiple of n. Interestingly 
this condition is the same as the condition for the occurrence of the nested peak structure 
in Sec. IIHI The comparison between the inset in Fig. ITTT a) for the large resonance peaks 
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and Fig. 0(b) for the nested peaks indeed verifies this relation between the two phenomena. 
Recalling that the inter-channel coupling is weak for iV c h < 10 and a < 10 x 10 -12 eVm, 
the correlation between the number of resonance peaks and N c ^ can also be explained easily 
since each channel separately gives rise to resonances. 

Lastly we examine effects of the nonideal spin injector /collector on a shorter SFET. 
FigurelT2lshows the conductance evolution for a SFET with relatively short L = 0.1 X 1.44 /im 
in a much wider window of a, (0 ~ 100) x 10~ 12 eVm. Again large resonance peaks appear, 
distorting the spin-orbit-coupling-induced modulation of the conductance considerably. We 
remark that the decrease of L clearly enhances the spacing of the resonances both in the G 
vs. a plot and in the G vs. Ep plot (not shown). However considering that the decrease 
of L also widens the minimum size of the a "window" for the observation of the spin-orbit- 
coupling-induced modulation of G and that the change of Ep simply "parallel" shifts the 
resonance positions in the G vs. a plot, reducing L does not appear to significantly facilitate 
the observation of the spin-orbit-coupling-induced modulation of G. 



B. Scattering effects 

Another practical problem ignored in the preceding sections is the scattering effects. 
Scattering by impurities, which are isotropic, short-ranged, and spin-independent, can be 
modelled by using the TB Hamiltonian of the impurity potential given by #w, ur i ty = 
Yl Sz Si=a Ylf=i Ui,j c tj Sz c i,j,sz> where Uij is a random variable whose value is uniformly dis- 
tributed in the range [— W, VF]. We investigate briefly the scattering effects in a weakly 
diffusive regime; W is chosen in such a way that the mean free path I is L/3. The nonideal 
spin injector /collector is again assumed and L = 1.4Afim, Ep = 0.398t are used. Figure IT3~1 
shows G as a function of a for various realizations of the random potential Uij. For all re- 
alizations, I is the same. Note that for certain realizations, the spin-orbit-coupling-induced 
modulation of G is clearly visible while for other realizations, the modulations is strongly 
distorted by the resonance peaks. This difference is due to the fact that the positions of the 
resonances depend on details of the spatial profile of the random potential U^j. A similar 
result is obtained for a nonballistic single-channel SFET j^. The results in Fig. 1131 also 
imply that the conductance of a nonballistic multi-channel SFET equipped with the non- 
ideal spin injector/collector may show considerable sample-to-sample fluctuations. Recalling 
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that there are more resonance peaks for larger iV c h, the sample-to-sample fluctuations are 
expected to be more significant for larger More systematic studies about nonballistic 
multi-channel SFETs with the nonideal injector /collector are necessary. 



V. CONCLUSION 

The issue of multi-channel effects were addressed. Via the numerical conductance calcula- 
tion using the tight-binding approximation, it has been demonstrated that the conductance 
modulation becomes less ideal with the increase of the channel number N c ^. For the ideal 
spin injector/collector, the conductance modulation ratio G max /Gmm as a function of iV^ 
(from a few to 90) has been examined for the RSO coupling parameter a ~ 5 x 10~ 12 
eVm and for the larger a ~ 50 x 10~ 12 eVm. It has been found that the decay of the 
modulation ratio with N c ^ is considerably faster for the larger a. Through the study of 
the energy dispersion relations and the spin configurations of eigenstates, it has been re- 
vealed that with the increase of iV c h and a, the spin configuration deviates further from 
the ideal ID configuration due to the formation of the phase-slip-like structure and the 
gradual drift of the spin angle, and that the energy dispersion deviates further from the 
ideal quadratic dispersion since the avoided crossing distorts the dispersion in a wider range 
of the longitudinal wavevectors. It turns out that those deviations are rather small when 
2m*aw/h 2 <C 27r while they become significant when 2m*aw/h 2 > 2ir. Thus this result pro- 
vides insights to the regime 2m*aw/h 2 > 2tt, where the inter-channel coupling is expected 
to be significant j^,^^. Effects of the electronic coherence in the multi-channel regime were 
also examined and we have found the nested peak structure, which was attributed to the 
coherent Fabry-Perot-like interference in the multi-channel regime. We have also verified 
that the magnetic field-induced peak splitting, reported previously for a single- channel 
SFET, persists in the multi-channel regime. For the nonideal spin injector /collector that 
consists of a conventional metallic ferromagnet-thin insulator-2DEG hybrid structure, it has 
been found that large resonance peaks due to the electron confinement by the insulators 
may strongly distort the spin-orbit-coupling-induced modulation of the conductance. When 
the 2DEG becomes weakly diffusive, the conductance modulation signal shows considerable 
sample-to-sample fluctuations. 
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APPENDIX A: IDEAL INJECTOR AND COLLECTOR 

As a test of the ideal spin injector, we calculate the electron transmission probability 
^Interface through the injector-2DEG interface for an electron incident from the injector side 
with the incident momentum k = (k x , k y ) and its spin pointing along the favored direction, x- 
axis. The model Hamiltonian of the injector-2DEG hybrid system reads f/hyb = ${x)H2D + 
$(— x)Hi n j, where *&(x) is the Heaviside step function. Since we are mainly interested in 
interface effects on the spin injection, the confinement potential V c in if 2 D and H [n j is ignored 
below. For a = 0, T interface is 1 since H 2 n and if in j are identical for the electron with its spin 
pointing along the x-axis. For nonzero a, Ti n t er face will deviate from one. The evaluation 
of T^terface is straightforward. Recalling that the eigenspin direction in the 2DEG differs 
from (+x)-direction and that the injected electron wave is thus a superposition of two plane 



n 

c 0.6 
.0 











to 


- - - T, 


to 


h 

— interface (=T^T 2 ) 

" fc)' 



30 60 90 30 60 90 30 60 90 
a(10" 12 pfm) a(10" 12 pfm) a(10" 12 fVm) 

FIG. 14: (Color online) Transmission probability Tj nte rface through the injector-2DEG interface as 
a function of the RSO coupling parameter a. The electron is injected from the injector side with 
its spin pointing along the (+x)-direction and with the incident angle 6 (0 for normal incidence), 
(a) 6 = 0, (b) = ir/6 and (c) 9 = it/ '4. T\ and T 2 (T\ + T 2 = T inter f ace ) are the transmission 
probabilities into the two involved eigenstates in the 2DEG. 
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waves with different spin directions, one needs to match properly the superposition with 
the electron state in the injector. The result is shown in Fig. El for the three incident 
angles 9 = 0, tt/6, and 7r/4, where 9 = sin - 1 (/c y /|k|). Note that regardless of 9, T interface 
remains close to 1 up to considerably large a, 100 x 10~ 12 eVm. Note also that T\ and T 2 
{T\ + T 2 = T inter f ace ) , transmission probabilities into each of the two plane wave states in the 
2DEG, remain close to their values cos (6/2 + vr/4) and sin 2 (0/2 + vr/4) in the limit a -»• 0, 
indicating that the spin of the injected electron points along the x-direction right after the 
injection. Figure El thus verifies the ideal spin injector. The ideal spin collector can be 
verified in a similar way. 

APPENDIX B: SYMMETRY ANALYSIS 



Symmetries often facilitate the analysis of a system and provide useful information 




Since if 2 D [Eq. is invariant under the translation along the x-axis, an energy eigenfunc- 
tion can be written as y) = e 1 [tpi(y)xi+^>i(y)xi)-i where ip^i(y) is the spin component 
along the (indirection, \-\ = ( I ) Pauli spin-up state and xi = ( ° ) spin-down state. 

i?2D also commutes with two other symmetry operators and Il 2: D S pi n (x, 7r), where is 
the time-reversal operator, 11^ is the mirror reflection operator with respect to the T/z-plane 
in the orbital space, and P sp i n (x, it) is the rotation operator with respect to the x-axis by 
the angle it in the spin space. To examine implications of these symmetry operators, it is 
convenient to use the representation £> spin (x, it) = exp(—i7rS x /h) = —ia x , and = a y K, 
where K is the complex-conjugation operator in the real space representation with the 
z-axis as the spin quantization axis. Application of 0IT c X> S pi n (x, it) on ^(x, y) results in 
e lkxX [ipt(y)xi — ipl(y)Xl\- Since [H 2 b, 0n x 'D S pi n (x, n)] = 0, the resulting state should also 
be an energy eigenfunction with the energy eigenvalue same as that for ty(x,y). Then by 
comparing V(x,y) with 0n x D spin (x, 7r)\I / (x, y), one finds that without loss of generality, 
ip](y) may be chosen to be real and ipi(y) may be chosen to be imaginary. Recalling the 
spinor representation y j f° r the spin with the polar angle 9 and the azimuthal angle 

<p, one finds that <fi for the energy eigenf unctions is ±7r/2 and that the spin direction of the 
energy eigenfunctions always lies within the ?/z-plane. 

When the confining potential is symmetric V c (y) = V c (w — y), which is the case for the 
hard wall confinement potential, if 2 D allows another symmetry operation zilj,P sp i n (y, 7r), 
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where U y is the mirror reflection operator with respect to the y = w/2 plane in the orbital 
space and T> sp - m (y, tt) is the rotation operator with respect to the y-axis by the angle n in 
the spin space. Since [iU y V spin (y , it)} 2 = I, possible eigenvalues of iU y V spin (y, tt) are +1 
and —1, and the energy eigenfunctions can be classified into even and odd parity states 
accordingly. For notational convenience, we choose <fi = +tt/2 for even parity states and 
= — it / 2 for odd parity states. At the special point y = w/2, U y plays no role and the 
parity becomes identical to the eigenvalue of iV sp i n (y, n) = a y . Thus at y = w/2, the spin 
of the even parity states points towards (+?/)-direction {6 = ir/2, = +vr/2) and the spin 
of odd parity states points towards (— y)-direction {6 = ir/2, <p = —tc/2). For y ^ w/2, 
the spin in general deviates from (indirection. But due to the parity constraint, the spin 
directions at y and w — y are correlated; 9(y) + 6{w — y) = tt. 

APPENDIX C: HARD WALL CONFINEMENT : EXACT RESULTS 

An exact energy eigenstate ^/(x, y) = e tkxX [tl'i(y)x j \+ipi(y)Xi\ °f H 2 r> [Eq. ((H)] with energy 
E can be expressed as a linear superposition of plane waves as follows, 



absence of the confinement potential [V c (y) — 0], r = (x, y), kj = (k x ,k y j), and k y /s are 
four solutions of E = h 2 (k 2 x + k 2 y>j )/2m* + (k 2 x + k 2 y j ) 1 / 2 (see Fig.©. When all four 

solutions are real (denoted by the blue solid line in Fig. 0|, one ma y t & ke k y> i = —k Vj3 and 
k V)2 = —k y ^. But in certain situations (denoted by the red dashed line in Fig. HJ), only two 
solutions are real and the other two remaining solutions are pure imaginary. We remark that 
in order to obtain a complete energy spectrum, one should consider such situations as well. 
In such situations, one may take k Vy i = —k y ^ for the two real solutions and k* 2 = —k y> 4 for 
the other two purely imaginary solutions. The procedure described below applies to both 
situations (blue solid line and red dashed line) denoted in Fig. HJ For the spinors Mi(kx) and 
■u 2 (k 2 ), we choose the representation, 



i 




(CI) 



where the plane wave (x, y) 



e^^Ujikj) is an energy eigenfunction of H 2 t> in the 




(C2) 
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and 

V2 - 



«*M = -i( • a L-, 2 I- ( C3 ) 



l«l 

For the spinors ^3^3) and u^k^), we choose their representation in such a way that the plane 
waves ^ kj (x,?/)'s satisfy the relations, ^ ka (x,y) = Hl y V spin (y, 7r)^ kl (x, y) and ^ k4 (x,y) = 
inj / r , spin (y, 7r)\I/ k2 (x, y), which is possible since iU y U sp i n (y, n) is a symmetry operation of 
the system. 

The four coefficients c^-'s in Eq. flCl|) are fixed by the hard- wall boundary conditions, 
ty(x, y = 0) = ^(a;, ?/ = w) = 0, which lead to a secular equation for a 4 x 4 matrix. Solving 
this equation is facilitated by the symmetry operation iU y U sp i n (y, n). Since iU y V sp i n (y , n) 
has only two eigenvalues, +1 and —1, the eigenstates can be classified into even and odd 
parity states. For an even parity state, one may choose C3 = c\ and C4 = C2, and for an 
odd parity state, one may choose c 3 = —c\ and c 4 = — c 2 . Thus the secular equation for the 
4x4 matrix is block-diagonalized into two 2x2 matrices, one for even parity states and 
the other for odd parity states. By solving the two 2x2 matrices separately, one obtains 
the eigenenergy E and the corresponding eigenstate *&(x, y) for even and odd parity states. 
Figure El shows the energy dispersion for various combinations of a and w. Figure |B1 shows 
the spin profiles of selected eigenstates. 

APPENDIX D: PERTURBATION THEORY 

Although the exact eigen wavefunctions and energy-dispersion relations are readily avail- 
able, perturbation theory calculations are still useful in order to gain insights into the a- 
dependence of eigen transport channels. Below we treat the RSO coupling as a perturbation 
and examine its effects on eigen transport channels. Since the longitudinal momentum k x 
is an exact quantum number of if 2D in Eq. one may work with the reduced Hamilto- 
nian H2T){k x ) — Hq + Hr that acts on the transverse part ip(y) of the original wavefunction 
y) = e xX ip(y). Here the unperturbed Hamiltonian H and the perturbation Hr are 
given by 



h 2 h 2 v 2 
2m* 2m* 



I, H R = a (vxjr - CTyk^ , (Dl) 



where I is the identity operator acting on spin. 
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In the zeroth order in a, each energy eigenvalue of H2n(k x ) is doubly degenerate with 



E^} = E$ = h 2 k 2 x /2m*+£ n , where E n = h 2 k 2 ny /2m* is the energy of the n-th (n=0,l,2,- ■ ■ ) 
excitation mode in the transverse motion and k n ^ y = (n + 1)tt/w. The corresponding eigen- 
functions are = 0n(y) • (x T + i Xl)A/2 and ipnl(y) = 4> n (y) • (x T -«xJ/\/2 respectively, 



,(o) 



are the spinors pointing along 



where n (y) = a/^T^ sinA^y, and \ ( i • 
the (±,z)-directions. 

Corrections to the eigenenergies can be obtained easily. The first order corrections to 
E^l (i =1,2) are given by the representation of Hr within the two-dimensional degenerate 



subspace: 



,(°) 



(0)| 



,(0)v 



,(o) 



(0)| 



,(0)v 



afc. 



(D2) 



Thus one obtains the first order energy correction \ = (—l) l ak x . Note that since the 2x2 
matrix of Eq. (jD2|) is diagonal, ip^ \ ari d Vn°2 are the proper zeroth order eigenfunctions in 
view of the degenerate perturbation theory The second order energy corrections are given 
by 



E 



(2) 
n,i 



V 



R- 



(D3) 



where the projection operator V n is defined by V n = 1 — J^ =1 |Vn°i)(V4°]l- By using the 
identity "P n a y iVw) = and the relation [H Q , o x \ = 0, Eq. ()D3|) can be written as 

2 



E 



which reduces further to 



(2) 
n.i 



a 
J? 



E 



(2) 



m*a 2 

"~2fi2~ : 



with the help of the identities p y = (m*/ih)[y, Ho] and ^ (f)*(y')(f>i(y) = 8(y — y' 
to the second order in a, the energy dispersion relation is given by 



Vk 2 x , _ 



-l) l ak x - 



* 2 
m a 



(D4) 

(D5) 
Thus up 

(D6) 



2m* ' " ' v _/ 2fr 2 

The result is shown in Fig. |5J Here we remark that although Eq. (jD6|) is derived for the 
hard-wall confinement potential, the first- and second-order corrections to the eigenenergy 
are independent of V c (y); its change alters only the zeroth order contribution S n . A similar 
result for a harmonic confinement has been reported previously 
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Next we examine corrections to energy eigenstates. According to the degenerate per- 
turbation theory, the first order corrections {ip^}) = ^iVn^) + Qn\ipnj), where Q n = 
E?=i and V n + Qn = l, are given by 

Wgi) = EElO (o) 1 , ) (D7) 



mj^nj=l 



E (o)_ E (o) 

n,i mj 



and 



(D8) 



where i = 2 for i = 1 and i = 1 for i = 2. Both Eqs. ()D7|) and (jD8|) can be consider- 
ably simplified. With the help of the identities V n ay\i/j^l} = 0, p y = (m* /ih)[y, H ], and 

E« tfiy^My) = &(y-y'), one finds 



ma 



(v-<v>»Wi°i>, 



(D9) 



Q»I^S> = o, (D10) 

where (y) n = J dy4>^(y)y(f) n (y) = w/2. Then up to the first order in a, the eigen wavefunc- 
tions are given by 

/w. . . m*ni(ii — \ /r,\ 

(Dll) 



^(y) = ^(y) + (-D^ m ' a(y fc 7 <y>w) tf}(y), 



We remark that although Eq. (jDll|) is derived for the hard wall confinement, it holds for 
other confinement potentials V c (y) as well when (y) n is evaluated for the new V c (y). The 



second order correction = 'Pnl'^ni) + Qn\4>n{) is given by 



(2)\ 



,(2)\ 



I n 



V 



It®, (D12) 



-^n.i n,i V n,i n,i 

With the help of similar identities, one finds 

2 



«2> 



p(o) _ p(0) , > , v 

•^n'.i ^V'.ii ,(0) \ \y)n",n' ty/n'.n 

WJ.i) 



n',n"j^n n,i n",i 



Qn\1>S/i) = 0, 



(D13) 



(D15) 
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where (y) n ,n' = (<f>n\y\<l>n') ■ We remark that these expressions for the second correction are 
again independent of the choice of V c (y). 

These expressions contain the information on the local spin direction. Since the Hamil- 
tonian if 2 D is invariant under the transformation QU x V spin (ic, it), the local spin direc- 
tion is always perpendicular to the demonstrated in Appendix [B] and it can 
be parameterized by a single angle 9 n ,i(y) within the yz plane; ip n ,i(y) = \i>n,i{y)\ • 
^cos ^i^-X] + i si 11 6n ' 1 ^ Xi 1 • The local spin direction is pointing — y, +z, +y, and — z axis 
when 8 n> i is — vr/2, 0, n/2, and ±7r, respectively. By using the perturbatively obtained wave- 
functions, the angle n ,i{y) can be expanded as a power series of a, 9 n>i = 9^J +9^J +9^J + • • • , 
where the zeroth and the first order contributions are 

€i(y)=-^(y-(y)n), (Die) 



= tfi(v) = -^(v -<*>»)■ (D17) 



and 



A similar result has been reported previously [25(. Note that the first order contributions 
change linearly with y, explaining the gradual drift in Fig. El Note also that despite the 
y-dependence of 9 ni the local spin directions of ip n ,i(y) and ip n ,2(y) remain antiparallel to 
each other \ijr n i(y)if) n ,2(y) = 0] at every position y up to the first order in a since \9 n ^{y) — 
0n,2(2/)| = ^ up to the first order in a. The second order contribution 9^\{y) can be also 
obtained (expressions not shown) from the perturbative expressions for ipn,i{y) and it turns 
out that when 9 n • (y)'s are taken into account, the local spin directions of ip n ,i{y) and ip n ,2(y) 
are not antiparallel to each other any more. Solid lines in Fig. El show 9 n ,i(y) calculated from 
the perturbative expressions for if) n> i(y), where it is clear that the deviation of |0 n ,i(2/) — 
0*1,2(2/)! from 7r becomes more evident for larger a. The second order contribution 9^J(y) is 
responsible for the phase-slip-like structure in Fig. El 



APPENDIX E: A NONIDEAL SPIN INJECTOR 



Here we consider the spin injection from a nonideal injector, which consists of a con- 
ventional metallic ferromagnet(F)-insulator(I)-2DEG heterostructure. The differences of 
the effective masses, the densities of states, and the Fermi wavelengths between the F 
and the 2DEG are taken into account and the electron transport within the F and the 
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2DEG is assumed to be ballistic. It has been reported 39] that the detailed structure of 
the spin dependent density of states within the insulator can influence the spin polariza- 
tion of the injected electron. In this paper, however, we ignore such details for simplic- 
ity and model the insulator by a 5-function potential US(x), where U denotes the tunnel 



contact parameter 
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41| . The model Hamiltonian of the F-I-2DEG hybrid system reads 



if^yb = $(x)H2t> + $(—x)Hf + U5(x), where Hp has the form Hp = p 2 /2m, 



electron ' 



■h a x + V . 



Here, ho represents the exchange interaction of the Stoner model in the F, and Vq denotes 
the difference between the bottoms of the energy bands for the F and 2DEG. Below we 
ignore the confinement potential V c and focus on the case of the normal incidence (9 = 0) 
from the F since the Fermi wave vector in the F is much larger than that of the 2DEG 29j 
and thus the normally incident electrons constitute a major portion of the injected elec- 
trons. We choose the value of the Fermi energy and the exchange interaction ho of the 
F to be 3.500 eV and 1.750 eV, respectively, and set the Fermi energy Ep of the 2DEG 
to E F = 0.398t = 0.103 eV. Thus V = (0.103 - 3.500)eV. Figure EH1 shows the spin po- 
larization P x of the injected currents as a function of the tunnel contact parameter U, 
where P x is a measure of the spin polarization along the (+x) -direction and is defined by 

( T inte = rLe ~ T inte = rfaee) / ( T inte = rLe + T inte = rfeee)- Here interface represents the electron transmis- 
sion probability through the F-I-2DEG structure for an incident electron with the spin s x . 
Note that while P x is only a few percent in the absence of U, it grows close to 30%. This 
demonstrates clearly that the spin- dependent tunneling over the insulator barrier is the main 
source of the spin polarization of the injected current while the bulk properties of the F are 
rather irrelevant factors for the spin polarization of the injected current. This observation, 



0.3 




U ( 1 0'^ J m ) 

FIG. 15: Spin polarization of injected currents, P x whose spin is along the (+x)-direction, as a 
function of the tunnel contact parameter U. 
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which is in agreement with Refs. |33j, |4l|, justifies the simplified description in Sec. IIV Al of 
the F-I-2DEG structure, where the differences between the two bulks (F and 2DEG) are 
ignored while the spin-dependent tunneling over the I is taken into account. 
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